0236 


; OT лға 


МАСА ТМ 3278 


а + 
-. 


з 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 


TECHNICAL NOTE 3278 


ATTENUATION IN A SHOCK TUBE DUE TO 
UNSTEADY -BOUNDARY -LAYER ACTION 
By Harold Mirels 


Lewis Flight Propulsion Laboratory 
Cleveland, Ohio 


Washington 
August 1956 


НЕА 
вт PL 


TECHN. 
МЇ її ` ғ 


! ! Ра 
ва wo жы, “тш лай 


— — ж a = 


5н29900 


ei 


-- -- -- 


| 


| 


р n mne 


TECH LIBRARY KAFB, ММ 


ІШІЛІШІІ 


TABLE OF CONTENTS 


Page 

INTRODUCTION. 52:34 9% de жож wo X^ 0 з 4 1 
ANALYSIS ©. м ; 5 
. Generation of Waves ^y Mass Sources : ; 5 
Generation of Waves by Uhsteady-Boundary-Layes "Action : 4 
ATTENUATION IN A SHOCK TUBE ...... 6 
Shock Attenuation Formulas .. oe ee әлі ш и йи 50% 6 
Limiting Solution Рог Weak Shocks $$ Else) ue de Us m пр 5,2 8 
Air-Air Shock Tube with T, = Т, = 5098........... 9 
RESULTS AND DISCUSSION ... gr de тй ме ысы Шо dee ДО 


Comparison with Experiments of Teferedce 6 E di ate а Ze 0 us lO 
Comparison with Theory of Reference 6 . . . . . . . . . . . . . 12 


CONCLUDING REMARKS ....................... 15 


APPENDIXES 
АБВОБ 5-ке ғ ЗӘР ЖР зш 2 ле Ate сш ЧЫЗ 


В - CHARACTERISTIC-LINE GEOMETRY AND SOME INTEGRALS ...... 18 


C - INTERACTION OF PRESSURE WAVES WITH INTERFACES . . . . . . . 21 
D - LAMINAR BOUNDARY LAYER BEHIND WAVE ............ 25 
E - TURBULENT BOUNDARY LAYER BEHIND WAVE ........... 29 
F - DEVELOPMENT OF ATTENUATION FORMULAS ............ 5 
G - IDEAL SHOCK-TUBE RELATIONS ................ 38 
Н - REYNOLDS NUMBERS AND TRANSITION. ............. 40 


H 
l 


GENERATION OF PRESSURE -WAVES BY WALL SHEAR AND 
HEAT ADDITION . . .. ; Я ғы е WO зи а wow ы. ES 


REFERENCES ........................... 46 


4059 


пе-1 


е 


NATIONAL ADVISORY СОММІТТЕЕ FOR AERONAUTICS 
TECHNICAL NOTE 3278 


ATTENUATION IN A SHOCK TUBE DUE TO UNSTEADY-BOUNDARY-LAYER ACTION 


By Harold Mirels 


SUMMARY 


A method is presented for obtaining the attenuation of a shock wave 
in a shock tube due to the unsteady boundary layer along the shock-tube 
walls. It is &ssumed that the boundary layer is thin relative to the 
tube diameter and induces one-dimensional longitudinal pressure waves 
whose strength is proportional to the vertical velocity at the edge of 
the boundary layer. The contributions of the various regions in a shock 
tube to shock attenuation are indicated. 


The method is shown to be in reasonably good agreement with exist- 
ing experimental data. 


INTRODUCTION 


A shock tube consists of a fluid at high pressure (region 4 of fig. 
1(а)) separated by a diaphragm from a fluid at low pressure (region 1). 
When the diaphragm bursts, a shock wave propagates into region 1 while 
an expansion wave propagates into region 4. А time-distance plot of 
these waves under ideal conditions is indicated in figure l(b). Regions 
2 and 3 have tbe same velocity and pressure but have different tempera- 
tures. The interface between regions 2 and 3 is referred to ав a con- 
tact surface. The analysis of the flow for perfect fluids is straight- 
forward (see, for example, ref. 1). Іп an actual shock tube, however, 
viscosity and heat conduction can not be ignored. These lead to a bound- 
ary layer along the walls of the shock tube as indicated in figure l(c). 
The boundary layer introduces nonuniformities into the shock tube. Ana- 
lytical studies of this boundary layer are presented in references 2 to 
6. One of the important consequences of the wall boundary layer is that 
it generates weak pressure waves which catch up with and attenuate the 
shock wave propagating into region 1. This attenuation has been studied 
experimentally and analytically in the work of references 1, 4, 5, and 
6, and is the subject of the present report. It is assumed that the 
boundary layer is thin relative to the shock-tube diameter. This is a 
practical restriction, since most shock tubes are designed so that the 
core of potential flow is relatively uniform in order to permit aero- 
dynamic tests. | 
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А few remarks concerning previous shock-wave-attenuation analyses 
are appropriate. In reference 4, the coordinate system is defined so 
that the shock wave is stationary. The flow between the shock wave and 
the contact surface is considered as a one-dimensional steady flow so 
that at each instant the mass flow through the shock wave equals the 
mass flow at the contact surface. If the mass flow at the contact sur- 
face is known at each instant, the corresponding shock strength can be 
found. The mass flow at the contact surface is determined from the lo- 
cal boundary-layer displacement thickness and free-stream conditions 
corresponding to an unattenuated shock. However, it can be shown that 
free-stream conditions do not remain constant at the contact surface 
(because of perturbations induced by the boundary layer). Moreover, the 
method of reference 4 does not take into account the existence of weak 
pressure waves between the contact surface and the shock wave. Since it 
is precisely these pressure waves which are responsible for shock atten- 
uation, the method of reference 4 can not be expected to give accurate 
quantitative results. It does, however, indicate some of the important 
parameters involved. Reference 1 uses the method of reference 4 in its 
studies of shock attenuation. At best, qualitative agreement with ex- 
periment is indicated. 


In reference 6, the flow perturbations due to the boundary layer 
are considered in detail. At each section of the shock tube, the veloc- 
ity and temperature variations associated with the boundary layer аге 
averaged across the tube to provide an equivalent one-dimensional flow. 
The wall shear and heat addition due to dissipation and heat transfer 
at the walls (all found from boundary-layer theory) are assumed to act 
on this equivalent one-dimensional flow, their action resulting in the 
generation of one-dimensional pressure waves propagating in both the 
upstream and downstream directions. By integrating along characteris- 
tic lines, the attenuation of the shock propagating into region 1 is 
then found. The theoretical trends appeared to be in good agreement 
with the particular experimental results reported therein. 


Reference 5 considers the flow in a coordinate system which is sta- 
tionary with respect to the shock. The unsteady nature of the flow be- 
tween the shock wave and the contact discontinuity is associated with 
the receding of the contact surface with respect to the shock wave. 
Weak pressure waves are assumed to be generated by the recession of the 
contact surface. The magnitude of these waves is obtained by a one- 
dimensional averaging procedure similar to that of reference 6. These 
waves overtake the shock and result in attenuation. The method of ref- 
erence 5 ignores the contribution of region 3 and can not be expected 
to yield good quantitative agreement with experiment. 


Of the previous reports on shock attenuation, reference 6 appears 
to give the best agreement with experiment. However, it can be shown 
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that the method of reference 6 does not give a completely valid repre- 
sentation of the wave phenomena induced by the boundary layer along the 
shock-tube wall. The deficiency of reference 6 is mainly associated 
with its use of wall shear in the determination of the perturbation 
pressure waves generated by the wall boundary layer. It is well recog- 
nized that the effect of a boundary layer on its external flow is di- 
rectly related to the vertical velocity &t the edge of the boundary 
layer. For example, reference 7, which is concerned with the Rayleigh 
(impulsive-plate) problem for a compressible fluid, shows that the bound- 
ary layer generates pressure waves in the external flow which are equi- 
valent to those which are produced if the wall moves normal to itself 
with & velocity equal to the vertical velocity &t the edge of the bound- 
ary layer. The nonzero pressure gradient over a flat plate moving at 
high speeds (because of the finite displacement thickness of the boundary 
layer) is an equivalent steady-flow phenomenon (e.g., ref. 8). Thus a 
proper way to find the waves generated by the unsteady wall boundary 
layer in a shock tube is to base the calculation on the vertical velocity 
at the edge of the boundary layer. Such an analysis is presented herein. 
The quantitative results for shock attenuation thus obtained would be 
expected to differ from those of reference 6. The comparison between 

the method of reference 6 and that of the present report is discussed 
more fully in the main body of the report and in the appendix titled 
GENERATION OF PRESSURE WAVES BY WALL SHEAR AND HEAT ADDITION. 


ANALYSIS 


One-dimensional flow with mass sources is treated, and the waves 
generated by these sources are derived. The waves generated by unsteady- 
boundary-layer action in a tube are then found, assuming the process can 
be considered as a one-dimensional unsteady flow. The application to 
the shock-tube problem is then indicated. 


Generation of Waves by Mass Sources 
Consider one-dimensional uniform flow in a tube of constant cross- 
sectional area. Assume that weak mass sources, uniformly distributed 
across each cross section, are ‚ and perturb the flow. Denote the 


net perturbation of a quantity (from the uniform-flow conditions) by A. 
(See appendix A for definition of symbols.) The equations of motion are 


(Ва Фа) eset] 


дар + р = +u са =m (Continuity) (1) 


Ap = a“Ao (Isentropy) 


— — —À M — M —— n ғғ 


———— — —— I — MÀ —— UM a ——— ет а --- aM aM 


eae aX а 


4. NACA TN 3278 


where m= m(x,t) is the rate of mass addition per unit cross-sectional 
area per unit x. The presence of the sources generates waves. Let the 
superscripts + and - indicate perturbations associated with waves 
moving іп the +x- and -x-directions, respectively. If &,t designates 
integration variables for x,t, the solution for the net perturbation at 
any point x,t may be expressed as 


Ap = Apt + Лр” (2a) 


= = (ар! - до”) (2) 
др = -5 Ap (ге) 


where 


The integrations are conducted along the characteristic lines 

т =t - (x-z)/(atu) and т-% - (Е-х)/(а-а) in the Ел plane. The 
upper limit on the integral for Ap” is ғо or -œ depending on 
whether М< 1 or M> 1, respectively. Equation (2b) incorporates 
the acoustic relations Apt = ра Ли" and Лр” = -pa Au. 


Generation of Waves by Unsteady-Boundary-Layer Action 


Consider a tube of uniform cross section to have flexible walls 
such that a small normal velocity у can be generated at the walls. 


This is equivalent to mass entering the tube at the rate pv dl, per 


unit x, where the integral is taken around the perimeter of the tube 
cross-sectional area. If the flow in the tube is considered as one 
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dimensional,+ the equivalent source strength is 


T 


„ру 
а 


(5) 


where а= 4A/ 1 is the hydraulic diameter, and v = қу Тһе ех- 


pression Рог Ар can then be written as 
Le 
1 - X 
а-ы 
x 
(4) 


x 
_ & 1. х= 5 
ар СЕЕ 


Similar expressions can be written Рог Au and Ap. Та the case of 
waves induced by boundary-layer action, the v in equation (4) refers 
to the normal velocity at the edge of the boundary layer. Note that a 
positive v results in compression waves, while a negative v results 
in expansion waves. 


lie the flow through the tube cannot be considered as one dimen- 
sional, it is necessary to consider each element of tube surface as an 
elemental wave source of strength proportional to the local value of v. 
The net wave strength at any point in the tube is found from an integra- 
tion over the entire tube surface. 


In the problem of attenuation in a shock tube, the shock wave is 
considered to be uniform, laterally, at each value of x, but to de- 
crease in strength with increases in x. Hence, a one-dimensional ana- 
lysis is permissable. Actually, small lateral variations of shock 
Strength exist, particularly near the walls, since the boundary layer 
immediately behind the shock exerts & tbree-dimensional effect on shock 
strength. These lateral variations are ignored herein. 


The shape of the shock wave was previously studied (ref. 9). It 
was assumed that region 2 was infinite in extent (i.e., region 5 was 
neglected), and the problem was considered as a steady-flow problem in 
э, coordinate system moving with the speed of the undisturbed shock. It 
was found that the shock wave assumes a parabolic shape, but that appre- 
ciable shock curvature is restricted to a region near the wall less than 
a boundary-layer thickness in extent. 
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ATTENUATION ТМ А SHOCK TUBE 


Equation (4) can be applied to find the attenuation in a shock 
tube. The details of the analysis are described in appendixes B to Е. 
The resulting formulas for shock attenuation are presented in the fol- 
lowing section. The limiting solution for weak shocks and a numerical 
solution for an air-air shock tube are also indicated. 


Shock Attenuation Formulas 


ЄСО? 


The flow in a shock tube is assumed to consist о? the ideal basic 
flow plus small perturbations due to the boundary layer. Ideal shock- 
tube flow relations are summarized in appendix G. The expansion wave 
of the ideal flow is assumed to have negligible thickness and to propa- 
gate into region 4 with velocity u = - a, ав indicated in figure 2. 
Let point d of figure 2 represent an arbitrary point on the shock-wave 
characteristic. The problem is to find the net pressure perturbation 
behind the shock (i.e., Аро, aJ." This requires an integration of equa- 


tion (4) along all the characteristic lines which contribute to Дро а. 


The major contribution to Ара, а comes from characteristic lines bd, 1 
be, and ab. Hence, these are the characteristic lines considered in the 
present analysis (as is indicated in appendix B). Integrating along 1 


these characteristic lines permits Аро а to be expressed ав 
2 


aod Ape а Үзг? 5% 


x Ес 
ШЕ С 
frg Yo = via Шм! vg ЧЕ + Ca А vs |х 
р 


Б 


-1 


1 - СЕ 202,с (5) 
Ар ,а 


Equation (5) is derived in appendix Е. The notation is defined in 
appendixes А and Е. The term СЕ Аро [es q represents the contribu- 
2 2 


tion to Дро а of all the characteristic lines other than bd, be, and 
2 
ab (appendix Е). 


ТР the boundary layer in region 2 is wholly laminar or wholly tur- 
bulent, the boundary-layer solution for və is known (appendixes D and 
E). Similarly, if the boundary layer in region 3 is wholly laminar or 
wholly turbulent, the solution for vz is known (appendixes D and Е). 
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For these special cases, the integrals of equation (5) can be readily 
evaluated. Thus, if the boundary layers in regions 2 and 3 are both 
laminar, equation (5) becomes (from eqs. (F6) and (F7) with 

nj = nz = 1/2) 


(Б ІШЕ ага = Ра, а _ -4y aLaF с Že : 
A - (u 1 - (us/22) 1/2 : 
| HM, и 


I 1 - e 
L+ — m 1/2 -1 
ея фы арт (5) | (6) 


8.52 l + Ма + age 


where Іс and Lz are found from appendix D. If the boundary layers 


in regions 2 and 3 are ке turbulent, equation (5) becomes (from eqs. 
(F6) апа (F7) with n, = = 1/5) 


а. V/9 га\/5 ара 5 дро а 72:50 YoDoF (s i 2 
(=) @ 9" pags didi 
Mp 


I 1 - (u 2] 
А. 4/5) -1 
Dr з2М2 Mp y/o 15 42. - сє (5) / ст) 


———————— (v 
azo 1 + Ма + 8.12 52 


where Го and Lz are found from appendix E. Similar expressions can 
be found for the ва cases 201/25 nz = 1/5, and Пр = 1/5, 
= 1/2. For other рои characteristics (such ав transition 


from laminar to turbulent flow in the middle of regions 2 or 3), it is 
necessary to integrate equation (5) with the corresponding ү distri- 
bution. Reliable criteria for determining the transition points in 
shock-tube boundary layers have not yet been established. A crude ten- 
tative method for estimating the transition points in regions 2 and 3 

is presented in appendix H. It is pointed out therein that the Reynolds 
numbers at point b, computed separately for regions 2 and 3, can be used 
as an index to determine whether the boundary layers in regions 2 and 3 
are primarily laminar or turbulent, respectively. The values of Ree b 


and Rez 4, as defined in appendix Н, can be obtained from figure 5 for 
2 
the case where the gas in region 1 is air at 520° В. 


- m+ ш ---- --“--- ee — — -.- н а-на —— I  —————=—— —————=——————ы=[—————=—————> 


8 МАСА ТМ 5278 


The present theory requires that the boundary layer be thin rela- 
tive to the shock-tube diameter. The boundary-layer thickness in region 
2 is defined by equations (D5) and (E3) for laminar and turbulent bound- 
ary layers, respectively. The laminar boundary-layer thickness бо їв 


taken to correspond to о/о, = 0.99, while the turbulent boundary-layer . 
thickness бо is taken to be the value obtained from a KÉrmán-Pohlhausen- 
type integral solution. Considering conditions along characteristic 
lines, the maximum boundary-layer thickness occurs at point b of region 

2. Values of бар апа бә ъ are plotted ав a function of М, in 
figure 4 for the case where the gas in region 1 is air at 520° R. 


It is previously noted, in connection with equation (4), that a 
positive v results in the generation of pressure waves, while a neg&- 
tive v results in the generation of expansion waves. From the 
boundary-layer theory of appendixes D and E it can be seen that, apart 
from dissipation and heat-transfer effects, vo is negative and vz is 
positive. Thus, the boundary layer in region 2 induces expansion waves 
(which attenuate the shock), while the boundary layer in region 3 induces 
compression waves (which accelerate the shock). Diss ipation and heat 
transfer modify these results. Dissipation tends to increase v in 
both regions 2 and 5. Heat transfer from the wall to the boundary-layer 
increases v. Heat transfer from the boundary layer to the wall de- 
creases- v. Та region 2 the heat transfer is from the boundary layer 
to the wall (ref. 2) and leads to larger negative values of vo and 


therefore more shock attenuation. Та region 3 the heat transfer is from 
the wall to the fluid for a weak expansion wave and from the fluid to 
the wall for a strong expansion wave (ref. 2). Thus, for a weak expan- 
sion wave, heat transfer in region 5 tends to generate compression waves 
(thus accelerating the shock). For a strong expansion wave, the heat 
transfer in region 3 tends to induce expansion waves (which attenuate 
the shock). The net effect of all these factors is to attenuate the 
shock wave. The relative magnitudes of the various terms in equations 
(6) and (7) are noted later in a numerical example. 


Limiting Solution for Weak Shocks 


For M, approaching 1, equations (6) and (7) take on the follow- 


Lamjnar case: 
1/ 
Us а224\ Ара ү с 
Up x ма = = ———————— |1 + (for Та ри = T4) 
a Vary + 1) Vo, 
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Turbulent case: 


Gy'Gj GE) ma. ың аур ба ey 
= pL À— Aue (for Т2 w = T5 r) 

(ез) (ap) 

9b 


The condition Tp w= Ty (eqs. (8а) and (9a)) corresponds to the case 


ere the chock- peim T dod perfect conductor (assuming that the 
fluid in region 1 is in thermal equilibrium with the wall). Note that 
taking T2, у equal to T4 means that the wall remains at its original 


temperature and that heat is transferred from the fluid (of region 2) 
to the wall. The condition T$ = To p (eqs. (8b) and (9b)) corres- 


' ponds to the case wherein the wall is a "perfect insulator or has & very 


low heat capacity. When the shock-tube wail is a metal, the assumption 


Te, у = Тт should give very accurate results, as discussed in reference 


2 (particularly since weak shocks are now being considered). Results 
for both То = T4 and T5 mm То „ 8ге given во аз to define the 
2 


upper and lower bounds of the heat-transfer effect on shock attenuation. 
The attenuation is greater, by a factor of approximately Yi; when there 
is heat transfer as compared with the insulated-wall case. (The factor 
is exactly Yi for the turbulent case and approximately үу for the 


laminar case, the latter depending on the value of оу). 


Air-Air Shock Tube with T, = T, = 520? В 


Equations (6) and (7) were evaluated numerically for an air-air 
shock tube with T4 = Ту = 5209 В. It wes assumed that о = 0.70 and 


Y = 1.4. Тһе results for the laminar and turbulent cases are summarized 
in figure 5. Both the insulated-wall case (т, Т Т; E Ts " = Ту, E 


and the case where the wall is a perfect conductor (т, nt is. v^ T. = Т) 


are noted therein. When the shock-tube wall is а metal, the assumption 
that the wall is a perfect conductor should give accurate results ex- 
cept possibly for very strong waves (e.g., ref. 2). 


The relative contribution to the net attenuation of the various 
terms in equations (6) and (7) is indicated in figure 6 for the 


= Т, case. In particular, the percent contribution 


———— — ae a te ee ee - 
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of the integrations along the characteristic lines bd, ab, and be and 

of the reflected wave at point c are indicated therein. (the reflected 
wave at point c represents the contribution of all the characteristic 
lines in figure 2 other than lines bd, ab, and bc, ав is mentioned in 
appendix F.) For weak shocks, the major ВВ to shock attenua- 
tion comes from the integration along line bd. With increasing М., the 
contribution of the characteristic line be increases gradually to a val- 
ue of about 30 percent at M = 6.0. The integration along line ab leads 


to campression waves (which tend to accelerate the shock) and therefore 
is negative in figure 6. Its value decreases to about -20 percent at | 
M, € 2.0 for the laminar case and to about -35 percent аб M, = 2.75 \ 
for the turbulent саве, and then increases with increases in M,. (The 
influence of characteristic line ab may be somewhat overestimated in ~ 
the present analysis because of the assumption of an expansion wave of 
zero thickness.) The contribution of the reflected wave at point : var- 
ies from а value of zero аб M, = 1 to about -5 percent at М„ = 

The neglect of the latter contribution appears reasonable for the al 
of M, considered herein. 


RESULTS AND DISCUSSION 


Equations (6) and (7) define the attenuation in a shock tube when 
the boundary layer is wholly laminar or turbulent, respectively. These 
equations are now compared with the experimental and theoretical results 
of reference 6. 


Comparison with Experiments of Reference 6 


Measurements of shock attenuation were obtained in the investiga- 
tion of reference 6 by using а high-pressure shock tube having а 1/6- 
by 1/8-foot rectangular cross section. Air at room temperature was used 
in regions 1 and 4. The air in region 1 was maintained at atmospheric 
pressure. Four sets of runs, corresponding to Ра = 4.061, 5.764, 


7.455, and 17.915 are reported. The results of these runs are summar- 
ized in figures 7(a) to (d), respectively. Some of the theoretical 
curves of reference 6 are included in these figures. The theoretical 
predictions of equations (6) and (7) are also indicated in figure 7. 
The latter were found from figure 5 with ру = 0117 and by assuming 


the wall to be a perfect conductor (since the shock tube had metal walls 
except for a pair of schlieren glass inserts). 


The Reynolds number at the contact surface for these tests is given 
directly by figure 5. The boundary-layer thickness at point b of region 
2 for the test conditions of reference 6 is plotted in figure 8. it is 
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assumed in figure 8 that the boundary layer in region 2 is wholly turbu- 
lent and that the wall is a perfect conductor. The boundary-layer thick- 
ness is presented in the form of a relative thickness 265 ъ/һ where h 
is the smallest shock-tube dimension normal to the flow and equals 1/8 
for the shock tube of reference 6. When 282 „b/b is small. compared to 


l, the assumption of a thin boundary layer is "valid, and the theory of 
the present report is applicable. It may be seen from figure 8 that 


25 p/n $ 0.20 for x < 12 for the test conditions of reference 6. 
Since 25 ъв tends to overestimate the effective boundary-layer thick- 
ness (because of the asymptotic manner in which the о approach 


free-stream conditions at the edge of the boundary layer) ,° figure 8 in- 
dicates that the theory of the present report is applicable for compar- 
ison with the experiments of reference 6. 


In figure 7(a), the first data point (x € 5!) is close Бо the the- 
oretical value for a wholly laminar boundary layer. The other points 
fall somewhat above the theoretical curve for a wholly turbulent bound- 
ary layer. The Reynolds number per foot at the contact surface for this 
case is Вер ъ/х = 0.4х106 (fig. 5). Using Re = 0.5x109 as & rough 


indication of the transition Reynolds number, &s discussed in appendix 
Н, it is reasonable to have the first data point near the laminar curve 
and the other points near the turbulent curve. 


The experimental results for M, = 1.442 аге given in figure 7(b). 
The Reynolds number at the contact surface is Reo ъ/х = 0.8х106. Тһе 


data agree quite well with the values for а wholly turbulent boundary 
layer. There is &.slight tendency for the points to lie above the the- 
oretical curve which might be attributed to the short length of laminar 
boundary layer directly behind the shock wave. 


The experimental results for М, = 1.518 (Ве, „/х = 1. 1x10) are 


given in figure 7(c) and are in excellent agreement with the turbulent 
boundary-layer theory of the present report. 


In figure 7(d), the experimental results (M, = 1.792, 


Rez, ъ/х = 2. 82105) fall somewhat below the theoretical predictions of 
the’ turbulent-boundary-layer theory of the present report. 


In general, the theory seems to agree reasonably well with experi- 
ment for the range of data considered in figure 7. Figure 7(d) indi- 
cates the poorest correlation and suggests that the present theory may 
underestimate the attenuation corresponding to large values of M,. 


2А properly computed displacement thickness would probably give a 
better estimate of the effective boundary-layer thickness. 


--——=——————— Bebo EB Ioue Pam an pt a e a t P git qp meo С 
+_———-+-+——-+——_—-—_—_———<—<—__——— раци, рар а а ap а р —— —- 


що ——— —— —————— A: - - 


12 МАСА ТМ 5278 


Comparison with Theory of Reference 6 


The present analysis assumes a relatively thin boundary layer, and 
the calculations are based on the vertical velocity at the edge of the 
boundary layer. The numerical value of v depends on terms which are 
related to the boundary-layer velocity profile, dissipation, and heat 
transfer. See, for example, equations (D3a) and (D3c). These terms 
correspond to the use of wall shear, dissipation, and heat transfer in 
reference 6. Dissipation and heat transfer play similar roles (quali- 
tatively) in the present analysis and in reference 6. However, it can 
be shown that the velocity-profile term in the equation for v has a 
sign opposite to that of the wall-shear term in reference 6 Гог charac~ 
teristic lines ab and be (e.g., appendix I). Therefore, the integra- 
tions along lines ab and bc, in reference 6, overestimate and under- 
estimate, respectively, the contributions of these characteristic lines 
to shock attenuation. If the same boundary-layer theory is used, the 
attenuation calculations of the present report and of reference 6 should 
be in quantitative agreement only in the limiting case M, + 1 (for 


which the contributions of characteristic lines ab and be are negligible) 
and for those values of M, where the errors due to characteristic 


lines ab and bc tend to compensate. 


Since different boundary-layer theories were used in reference 6 
and herein, it is not directly possible to separate discrepancies be- 
tween the two methods due to the respective boundary-layer theories from 
the discrepancies due to basing the attenuation calculations on wall- 
shear, dissipation, and heat-transfer terms rather than on v. An 
estimate of the latter discrepancy can be obtained by reversing the sign 
of the first term in the equations for То (eq. (D3c)) and I4 (eq. 
(D4c)) for the integrations of v along characteristic lines be and ab. 
Such a procedure shows that the theory of the present report and that 
of reference 6 agree at М. = 1. With increasing M,, reference 6 first 
overestimates the attenuation because of the increasing importance of 
characteristic line ab. At М. = 0(1.5), reference 6 overestimates at- 
tenuation by about 10 to 15 percent. With further increases in M,, the 
errors in characteristic line be become important and tend to compensate 
for the errors in line ab so that the discrepancy between the two methods 
decreases. At M, = 0(2), the two methods are again in approximate 
agreement. As M, increases further, line &b becomes relatively 


less important compared with bc, and the theory of reference 6 underes- 
timates the attenuation. At M, = 6, the method of reference 6 appears 


to underestimate the attenuation by about 50 percent. The above fig- 
ures are only approximate, because of the manner in which they were ob- 
tained, but they indicate the proper trend with M,. 
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From figure 5(b) it can be seen that the attenuation theory of ref- 
erence 6 (for turbulent boundary layers) is in good agreement with that 
of the present report for 1<M,< 2. .The agreement is somewhat better 
than that to be expected from the discussion of the previous paragraph 
and is due, in part, to the use of ап incompressible-boundary-layer 
theory in reference 6. Hence, either theory could be used to correlate 
the experimental data of figure 7. For М, > г, reference 6 consider- 


ably underestimates the attenuation. The laminar-boundary-layer theory 
of reference 6 differs by a factor of about 2 from that of the present 
report; therefore, there is a large discrepancy between the theoretical 
laminar-boundary-layer curves in figure 7(c). 


CONCLUDING REMARKS 


A method is presented for computing the attenuation of a shock wave 
due to unsteady-boundary-layer action. The various assumptions involved 
in the analysis are summarized in this section, since these define the 
limitations of the method and suggest possible fields for improvements: 


1. Small perturbations. The equations of motion were linearized 
assuming the potential flow external to the wall boundary layer under- 
goes only small perturbations. For long shock tubes with large amounts 
of attenuation, it might be advisable to employ a characteristic method. 


2. Thin boundary layer. The assumption of a thin boundary layer 
(relative to shock-tube diameter) is consistent with the assumption of 
small perturbations of the potential flow. ТЕ the viscous effects span 
the entire tube cross section (i.e., long tubes), it may be advisable to 
base the shock-attenuation theory on wall shear. 


3. One-dimensional longitudinal waves. It is previously pointed 
out that the perturbations in а shock tube should be computed by assum- 
ing each element of wall surface area to be an elemental acoustic source 
of strength proportional to v. This gives rise to a complex wave pat- 
tern involving both longitudinal and transverse waves. For the purposes 
of the present analysis it was assumed that the longitudinal waves are 
of primary interest and that these can be computed on the basis of a 
simplified one-dimensional theory. This assumption is accurate when an 
Observer is relatively far fram the sources (since the details of the 
source distribution around a perimeter then become relatively unimpor- 
tant) but introduces errors when the observer is near the sources. It 
would seem that the assumption of one-dimensional longitudinal waves is 
accurate for weak shocks but introduces errors for stronger shocks, par- 
ticularly when the flow relative to the wall is supersonic. 
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4. Boundary-layer theory. The boundary-layer theory of reference 
2 was used. It may be assumed that the laminar-boundary-layer solution 
is reliable except possibly for very strong shock waves. Та the latter 
case, very large temperature gradients exist normal to the wall, and it 
may be advisable to choose a different reference temperature from that 
used herein. Also, for strong shock waves, dissociation might occur 
which would also require changes in the laminar-boundary-layer theory. 
The turbulent-boundary-layer theory of reference 2 requires experimental 
verification for even the weak~shock case. However, the good agreement 
between the attenuation calculation based on the turbulent boundary layer 
of reference 2 and the experiments of reference 6 suggests that the 
turbulent-boundary-layer theory of reference 2 gives reasonable results, 
at least for the weg&k-shock case. 


BGO? 


o. Expansion wave of zero thickness. In order to simplify the prob- 
lem of determining the boundary layer behind an expansion wave, it was 
assumed in reference 2 that the expansion wave is of negligible thick- 
ness (i.e., "expansion shock"). This assumption is valid for weak ех- 
pansion waves but is in error for strong expansion waves. The boundary- 
layer solution for region 3 might be improved for the strong-wave case 
by considering the finite thickness of the expansion wave. However, the 4 
contribution to shock attenuation of region 3 becomes small for the 
strong-wave case, so that an improved boundary-layer solution for region 
3 may not significantly affect the attenuation calculation. x 


In addition to the preceding discussion, the following extensions 
of the present report might be pursued: 


1. The details for obtaining the perturbations at a fixed point in 
the shock tube, as opposed to finding the shock attenuation as is done 
herein, might be treated. This would be useful for further correlating 
the theory with experiment and for determining conditions at ап аего- 
dynamic model (when the shock tube is used as a wind tunnel). 


2. Formulas equivalent to equations (6) and (7) might be derived 
for the case where the boundary layer is partly laminar and partly 
turbulent. 


5. The boundary-layer theory as outlined in appendixes D and E 
might be simplified, particularly for the turbulent case, so as to give 
reasonably accurate results without too tedious a boundary-layer 
calculation. 


Lewis Flight Propulsion Laboratory 
National Advisory Committee for Aeronautics 
Cleveland, Ohio, April 24, 1956 К 
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APPENDIX A 


SYMBOLS 


The following symbols are used in this report: 


--т — a - -am e auam = чат y що 


cross-sectional area of tube 

speed of sound 

eq. (C3) 

specific heat at constant pressure | 
eq. (C3) 

hydraulic diameter, 4А/1 

eq. (C7) 

eq. (C7) 

smallest shock-tube dimension normal to flow 

appendixes D and Е 

perimeter of shock-tube cross section 

Mach number of flow relative to wall 

shock Mach number relative to wall 

Mach number of flow in regions 2 and 5, relative to wall 
rate of mass addition per ait ares es unit x 
&ppendixes D апа E 

pressure, 1Ъ/ва ft 

Reynolds number (appendix Н) 

Reynolds number &t point b of region 2, 


&1X uo 2 Us | 
У. У] 2. a, 1 + Мо - a4) "в > appendix Н 


e wot на Tu ——— =©Ө'_—--—-————-—-—-——-_+_——-—_—-—-+-- —--—-—-—-——-—-—-— —-—-+:=—-—-—-—-—##. . —.—-—.——-——.—<—<—.—=——<—-— ami s Ra e — ------. 
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T temeraria; OR 

T. temperature of insulated wall 

% time 

u velocity of flow relative to wall 

us velocity of shock wave relative to wall 

Чо Ug velocity in regions 2 and 3 relative to wall 

ү vertical velocity (positive when directed into tube) at 

edge of boundary layer 

x longitudinal distance 

Y ratio of specific heats 

^ perturbation quantity (Ap = perturbation of р, etc.) 
` $5 laminar-boundary-layer thickness in region 2 (eq. (D5)) 
5, turbulent-boundary-layer thickness in region 2 (eq. (Е5)) 
u coefficient of viscosity 

v kinematic viscosity 

Е integration variable representing x 

р . mass density 

б Prandtl number 

д integration variable representing t 

T. wall shear stress, positive in 4x-direction 
Subscripts: 

1,2,5,4 regions of shock tube (fig. 1) 


@,b,c,d,e points on characteristic lines (fig. 2) 


m evaluated at mean temperature of boundary layer 
(appendixes D and E 


W evaluated at wall 


сепъ 


4059 


CC=3 


NACA TN 3278 17 


Superscripts: 

+ associated with wave moving downstream (+x-direction) 
Е: associated with wave moving upstream (-x-direction) 
Special notation: 


ато = &z/855 P12 = рү/Рә, etc. (Two successive integer subscripts, 
| not separated by а comma, represent а ratio.) 


— аз = a тт ———— < ————————— 
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CHARACTERISTIC-LINE GEOMETRY AND SOME INTEGRALS 


The equations of the characteristic lines considered in the present 
report are summarized. Some important integrals are indicated. It is 
assumed that the characteristic lines are straight (corresponding to the 
ideal-flow case) and that the expansion fan can be considered as an 
"expansion shock" (i.e., expansion wave of zero thickness) moving with 


velocity а, into region 4 (following sketch). The point а with coor- 


dinates x,t is assumed to be on the shock-wave characteristic. 


x,6 
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19 
The equations of the characteristic lines are: 
Line bd 
Е =x - (а + uj)(* -т) 
et oo 
82 + uo 
TOUR 
от -Е uM TX, sl (x - &) 
Line &b 
5 = ё - (az + uz)(%, - +) 
(B1) 
1 + Ма +a 
5 45 
ат +E = C PRA" (E - Ea) 
Line be 


E = Ey + (ар - ug) (тъ -т) 


&- Е 
Wo = мд 


1 -М 
ит -& Е, (E, -Е) 


LLL M 


20 МАСА ТМ 5278 


In terms of x,t, the coordinates of points а, b, and c are 


Point а 


қ B EIE 


e 
и 


: Е + Мо - быа) 


Point b 
ug 
Ey = у Ё + Мо - (а/а) | a, 


of 
! 


5 [1 + ме - (а/ш) (вг) 
Point с 


Е + М - Е] 


Sc TUS l - Mo + (u,/as) 


E. Ц (о/а) 
l - № + (u,/ag) 


vot 


Some important line integrals are 


x 


1 -no l-no Ug > 4" 
ТИ (ut - &) а& = (1-ug/u,) с + № - z) 1 - по 


Sp 
Ес по E ч 1-02 -no 1-02 
1 с | 5 x 
a; | (ат - Е) dë = =: - =) с + Mp - =) 1 по 
l-n 
u a, 5 
1 " 3 5 ( + Me- 4 4 di. 
TM; / (a+ E) 45 = (1495/84) 1 + Ма + ал» I - nz 
Ea 
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APPENDIX С 


INTERACTION ОЕ PRESSURE WAVES WITH INTERFACES 
Weak pressure waves are assumed to overtake а contact surface or а 
Shock wave. The strength of the reflected and transmitted waves is 
found. 


Contact Surface 


Consider the contact surface which separates region 3 from region 
2 (following sketch). At a certain instant, known incident waves AD y 
2 


+ = oe surface 


and Аро ы intersect the surface. The problem is to find the final 
2 
waves Арз b and Ape р. From isentropic flow relations, 


Ар? = pa Aut 
(c1) 
Ар” = -pa Ли” 
The boundary conditions across the contact surface are 
+ = 
А05 р + Ало = Aus р + Аш» ү T€ 
CZ 


I nde did e i PD——O————— E M — 
— SE ——=ыы—.-—==—=ыы—ы—=ы= == = +———— ——_—--————_—-+——— .——_—-—-—————.. sage 
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From equations (Cl) and (C2), the expressions for Ару p end Ap; 
are 4 


Е + 
Ары = Pe + Пару y 
(сз) 


- E 
Арз = T3223. р 7 5955 


where 


- -6GOP 


C= (туар - 1)/(f 32223 + 1) 


Ш 


2/28: +1) 


р 


The quantities С and D are reflection and transmission coefficients, 
respectively. 
Shock Wave 
Consider a shock wave moving with velocity us into & stationary " 


fluid (region 1). The region behind the shock is designated as region 
4. Conditions in regions 1 and 2 can be found from shock-wave theory. 


In particular, with Yi = Yo» 


ФМ - (vi - 1) 
Po. = TES! 


2. [мт 
аа Yi + lis Mg 
At & certain instant, & known wave Аро а intersects the shock (follow- 
2 


ing sketch). 
— discontinuity 


/_—-Shock wave 


(C4) 


қаныға 
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25 
Тре problem is to determine the net perturbation of the shock wave. 
From equations (C4), 
sr 
Apo а ЖҮ 
Po Yi +1 le 8,8 
(C5) 
^u» q 2 1 
&4 "Y, +1 GILT 
1 5 


The net pressure and velocity perturbations in region 2 are then related 
by 
5 


М5 
Ар, a = 2018) тм Aue а (ce) 


With Ap a? др, а + Ару, а and Аар, a= Аид) а + Au, a using equations 
(C1) gives 
- na + 
Ар а ^ ЕАР» а 
А (ст) 
“Pea ҒАР а 


3 


where 


5 MS 
артат 5 - 1 212842 в + 1 
12 214. Me 12 2 j + м 
5 


гл) 
| 


M. 


доза + 1 
2P10 ЖҰРТЫ 


Е = |4р12842 


1 + Me 


The quantity E is а reflection coefficient. The resulting perturba- 
tion of the shock Mach number is 


+ 
үз + 1D Ар 
за АҒ c yd (св) 
Б Yi ма P2 | 


———— --- ж. ——Ъ———— ст а ei ee ep ep ee аа аа 
SE Re E n rne 


- —— — 
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For М. =1, сова = 0. For M, infinitely large, 


TE 
_ TEE nod 
A2 d _ АР? а 


———————, which becomes ---2- = - 0.14 for Yi = 1.4, 


102 а пра NS Ар? a 


indicating that pn reflected wave is relatively small and is opposite 
in sign to the incident wave. 


COP 
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APPENDIX D 


LAMINAR BOUNDARY LAYER BEHIND WAVE 


The laminar boundary layer behind a shock or thin expansion wave 
is analyzed in reference 2. Some of the results are summarized herein. 
It is assumed that the wall temperature behind the wave is constant and 
that о and ср are independent of temperature. The fluid properties 


и and k are referenced to a mean temperature as discussed in appendix 
C of reference 2. The notation of reference 2 is compared with the no- 

tation of the present report in the following equations. The left-hand 

side of each equation represents the notation of reference 2, while the 

right-hand side is in the notation of the present report. 


Region 2 
E Uw 1. 
Us 7 Ug wu. I- АТ 
‚ Че 1 - учуш 
Че = U, - и, 
(D1) 
Uy ЖЕ по ав 
MB = Мо Ue “1 - а/а.) 
х= Uu - Е 
Region 3 i: 
Uy „=. 
Ч; = 84 ü l + (uz/a,) 
Че = а, Ж uz - | (D2) 
Ve = YZ Че 71+ 8, 


X = ат + Е 


Vertical Velocity at Edge of Boundary Layer 


The expression for the vertical velocity at the edge of the bound- 
ary layer can now be written as follows: 


Region 2 


тей b 
Е 03/ üs f v2 г (рза) 
ата тт гел СЫ: 


У—=—ы———-- —— A ge a re д 


(ет) - 52) 
Lim <<" “энд = 1.134 А (рза) 
|-> u2 “s 2.022 - (а/а) | 
-0.045u9/u, 
| 1-(u2/u,) 
со) (D3e) 
-0.5040.47 (us /u,.) 
1 - ( fu ) 
279 (зе). 
0.37-0.39(u,/u,) 
l - us 
r4(0) = (og) (255) 
To r У2 - 1 j 
ez „ша „в-ф вт) вола 2.2) (ова) 
22 ш 2 2 2 


The above quantities can be evaluated if the wall surface temperature 
То,м is known. The value of Tow depends on the heat capacity and 


conduction properties of the wall anü'the wall thickness. А method for 
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evaluating To м is presented іп reference 2. From the discussion іп 


reference 2, it may be concluded that in most cases То, w BSPDPproximately 
equals I4 when the wall is & metal. 


Region 3 
Б 12 Е l-nz vz nz 

та те BA (x) (D4a) 
D where 
© | 

ns = 1/2 | (D4b ) 
1 НЗ, що а + ee с 1 
2” E = —2- Мама (1 rs dn + 
т приза О 
ға 5 
о 
; 1 + (us/a4) (Т, 13 r j 
3! - ——- 55 91 (D4c) 
су u /а Tz Tz . | 
0.37+0.50(uz/a, ) 
1+(uz/a,) 
rz(0) = (oz) (D4g) 
ә т T = 1 г 
5, г 9 
Е = 1 + ——— Ме r.(0) (раһ) 


-- 4 ee T —RH- HU oe — ees nem nnnc der» ere а е хш —— —— — — — 
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T T Tu 
СЗ 0524 1 4 0.22[ ez _ (041) 
Ts У Ts 


Equations (D4) can be evaluated if Tz , is known. If the wall is in- 
itially at temperature T,, then Ts, y = Та ів generally a good esti- 
mate for Ts, y» 25 is discussed in reference 2. 


Boundary-Layer Thickness in Region 2 


In the body of the report it is assumed that the boundary layer is 
thin relative to the shock-tube diameter. It is therefore of interest 
to present an expression defining boundary-layer thickness. Region 2 is 
of particular interest since the boundary layer in region 2 is generally 
larger than that in region 3.. Let 6) represent the киык ие 


layer thickness in region 2, defined во as to correspond to — = 0.99. 
From equation (18) of reference Э, ul 


vo(u,t - Е) 2.89-/1 - цо ас е Ас. „12 - 1 2 . 
fase e edP eene 
1.415 - ufu. 
О 


=) вә dm | i (D5) 
1/0 | 


(Еа. бе) uses а mean reference temperature and an interpolation formula 
for 15 (defined in ref. 3) and takes the upper limit on the integrals 


to be » rather than 15.) To evaluate Бо along the characteristic 


1 + Me xe) TN 


line bà of figure 2, take. от - Е equal to DOES 


The value of бо at point b is found by replacing uT- Е by 


1 +M - u,/a u 
ы BE + M. 2) - x) x. 
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APPENDIX E 


TURBULENT BOUNDARY LAYER BEHIND WAVE 


The turbulent boundary layer behind a shock or thin expansion wave 


is also studied in reference 2 for the case of constant wall surface 


temperature. The correspondence between the notation of reference 2 
and that of the present report was previously noted in equations (D1) 
and (D2). The turbulent-boundary-layer solution of reference 2 (relat- 
ing to v) is now summarized. It is assumed that the wall surface tem- 
perature is essentially constant. ` 


Vertical Velocity at Edge of Boundary Layer 


Region 2 
| "ай, UT, у; NP 
уе = - № а Ca uin 
where 
по = 1/5 | (ЕТЬ) 


58 Cuz/us) - 1| Бао) 4/5 
5. (ug/u.) г 65 [(u2/u,) g 1] 


52 (ufu) -1 за. T ча 22 
ы буйр | (еол) (0-48) (өш 


| T2 Up Ug 
92 (ug/ug- 1) _ "Taw | 2567 (2 + ug J 27 * uy 12,8 ші 
бо o nm 1 - (up/u, ) 


ы (zie) 


Lo = 0.0460 
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1 
zN az 
i2,N == 2 (Elg) 
О 1 + 022 - Coz 
T n 
2.0 
Ъ2 = m. -1 (E1h) 
эм 
Р.Т To 
со =| ——-1 (Eli) 
с ( To Хе) 
To т у2 - + 1/5 
Filius Ji UE 
m s 1+7 ма (ог) (Е13) 
To m To м Т2 г 
A = 0.5 | 2 + 1 + 0.22 | -т- - (Elk) 
+O 2 | 2 


The conditions under which To y essentially equals Ту can be estab- 
lished by the methods of reference 2. The assumption Tow = Ту appears 
reasonable TOR most cases. 

The integral 12 А is tabulated in reference 2 Рог М = 6, 7, and 
8 and various values of bo and Co. The integral can be evaluated 
analytically from the following expressions: 


1 
ROC REX 21 22 зі І ) 
2,5 2 [ 2N B2 
ү + 2с, 


- са, | (E17) 


where 
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Под” y КЕ : Э (5) | (кір) 
‚ Ш 


Alternate ааа for (15 ww) Bo and (12 wr» are 
(т) = — Nim! 1 Ww 
га Ва  1.- Bo (Nel)! \ 1- Bo 
m=O 


E 1 па м 2! сеза) 
C CNN-a) |27 Ва) арма 


2 2 
‚ (+3) (+2) (1-8, ) 


(I5 yds = GOAT [ x (N Р 1) 2, Gays] (=) 


қ N+1) №-1. | 
“тайга | "(uu M Hs | 


(Elr) 

Region 3 

2 l-n n 
Уз, = L. 95/4 ` eom - (Е2а,) 
l + (м. /ал) алт + Е 
where 
n, = 1/5 | | (E2b ) 
ы 4/5 | 
8i |1 + 5 

Lz = 0.0460 вр + (5/54) ] ӛз) | (Е2е) 


5.(а./8,) 95 Əz |1 + (uz/a,4)| 


x | 
55 (1 + (us/a,) | ад uz Tz Uz 
| a 400554 a 
б де: 4 Ота Од M 


5 4 


= ————- ———- =—_————> m, арн mo мз м op gh t 
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ар + (us/e,)] _ т(Фз/Т» w) _ | _ 5 I - 55 І (Еге) 
быы) 1+ Cus), 15,6 MEO MN 
1/4 ` 5/4 
H3 m / 15 / 
o4 = (521 3 (E2f) 
Hz 5,ш 
т М 
2 dz 
AE —— - | (E2g) 
О 1 + 5-2 Б с-2 
bs = (Tz fms w) -1 | (E2h) 
ccm LES н || (25/75) е 
т Y - 1 ; 1/5 
5,г 5 
г. 
m m М т . 
x - 0.5( RR ЕА ) Е e = = ) (в2к) 
5 5 | B 


Тһе conditions under which Tz „ essentially equals T, can be estab- 
2 
lished by the methods of reference 2. The assumption Tz = T, appears 
2 
reasonable for most cases. The integral ls, y cen be evaluated from 


the tables of reference 2 or equations (E12) to (Еіг) "- the sub- 
script 2 replaced by the subscript 5 therein). 


Boundary-Layer Thickness in Region 2 


The symbol б represents the turbulent -boundary-layer thickness 
as obtained from an integral (K&rmán-Pohlhausen) type of boundary-layer 
solution. From reference 2, it сап be shown that 


СГ» 
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ғғ 


52 = 0.057420. 


(ЕЗ) 


The value of б along characteristic line bd or at point b may be 
obtained by the substitutions indicated after equation (D5). 


ee ea A trm IES mm жы MURR — ------ 908 7 
- -—- NNNM Annon Lem n аат. ц ци en рч она 
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APPENDIX Е 


DEVELOPMENT OF ATTENUATION FORMULAS 


The attenuation of the shock wave in a shock tube is now derived. 
The flow is considered to consist of the ideal shock-tube flow plus 
small perturbations due to the boundary layer. The ideal flow is uni- 
form in regions 1, 2, 3, and 4, and is denoted by the appropriate sub- 
scripts (Чуб; Py» Ро» Ра; etc.). Perturbations at any point are de- 
noted by A and an additional subscript indicating the point. Thus, 
the perturbations at point b of figure 2 evaluated on the right-hand 
side of the entropy discontinuity are designated by АР t Р Aue ъ ; апа 
so forth. The expansion wave of the ideal flow is assumed to have neg- 
ligible thickness and to propagate into region 4 with the velocity 
u=-a, as indicated in figure 2. (The same assumption was used in 
the attenuation study of ref. 6.) 


6407 


Let point а of figure 2 represent an arbitrary point on the shock 
wave. The problem is to find the net pressure perturbation behind the 
shock (i.e., Ap, а}: This requires an integration of equation (4) along 

2 


all characteristic lines which contribute to Ар, 4. Because of the 


entropy discontinuities, there are an infinite number of line segments 
along which the integration must be conducted (fig. 2). However, the 

major contributions to APs а can be shown to come from segments ab, 

cb, and bd. ы 


ТЕ Ар; а represents the incident pressure wave at point а, then 
equation (c7) gives 


re | 
Ар» аы. | (F1) 


But, from equation (4), 


| x 
+ + 
Аға = Poy ват + [o үг(Е,9)аЕ (F2) 
ur ща. ж | | 


where the integration is conducted along the line bd. From equation (C3) P 


+ 
+ D^p F3 


+ = СЛ 
“Ро СР, 
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Again, from equation (4), 


| PY 2P2 ж 
A95 ъ = Ар, с + ш | У2(Е,т)аЕ 


2г ар 

+ + 5 

б,” P3,a * a dl + М.) au та (Е s Jae 
ба 


It can be shown that E = О in the present case. Moreover, 
Ро = Ра: The expression zi Ара, а сап then be expressed as 


x : 
Paca NN NONE ТРИЕ ый aE А 
Е 2 Ро 1. 2 1-М 2 E ) 
5p En 
hs -1 
1 CE Фос 
2р2 а 


In appendixes D and Е it is shown that for wholly laminar ог wholly 


turbulent boundary layers 


92/84 sns Уз 05 
Vz = lz ee 
9 1 + ад (= + z) 


59 


(F4) 


(5) 


(F5) 


where L and n are independent of т and E. The values of L and 
n depend on whether the boundary layer behind a given wave is laminar 


сары ч— ee e eee ee ee ee eee ma 
———— nidi cl pe ай 


х 
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or turbulent. Substituting equations (F5) into equation (5) and using 
equations (B3) yield 


2-5п n 1- | -2по 
(%) (=) (2 Зада, Pra аа) | к 
чо Vo x Po 1-02 | 1-(а ДА j 2 х 
1 sf 2 
c 


l - CE —2— | (F6) 


The corresponding change in shock Mach number can be found from equation 
(C5). Appropriate values for Г and n are given in appendixes D 
апа E. 


If рә = п5 and СЕ Аро, c/ Ap. a 18 assumed not 5 vary with x, 
then equation (F6) indicates that Аро, с/Авг, а = (Е г/%) 780 апа 


l-n 
Аро : 2 
Ро d T" 
р 


Equation (F7) indicates СЕ Ape c/ A5. а Бо be independent of x and 


is therefore consistent. with the original assumption to this effect. 
Substitution of equation (F7) into equation (F6) gives accurate results 
even when по £ nz; since the major contribution to Ap» а comes from 


the first term on the right-hand side of equation (F6). "The term 
СЕ Ара [Po а represents the contribution to Аро а of all the char- 
2 2 2 


acteristic lines of figure 2 other than lines ab, bc, and bd, and is ге- 
ferred to as the contribution of the reflected wave at point с in fig- 
ure 6. It can be seen from figure 6 that this term is small and can 
probably be neglected in most cases. 


в 
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If it is necessary to consider the-boundary layer behind the shock 
{or expansion) wave as partly laminar and partly turbulent, then equa- 
tion (5) must be integrated accordingly. 
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APPENDIX G 


IDEAL SHOCK-TUBE RELATIONS 


. .Shock-tube relations, assuming ideal flow, are presented herein 
for convenience. The formulas were obtained from reference 1. The 
notation ру, = р/р, ‚ R35 = а „/а.) and so forth, is again used. 


Define 
а = (y + 1)/(r - 1) 
В = (ү - 1)/2 
Then 
1 
рул = Б. (Gla) 
ifr, 
Pz, = (PyyPoq) (алъ) 
1 + 91221 
Pan a FP = 
2В 
Ту, = азд = (P14P21) - (ала) 
Po, (a, + Ро.) 
_ 2 Рәт: * P21 
Тә] = 821 = I 1 + a P21 | (Gle) 
Ug Poy ~ 1 | 
a, Г. 1/2 (ва) 
Yal Вуболрру + 2 
us x d 
x ec eed (Glg) 


LONE 


муу 
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чә P23 74 


ет я ra [Boos b Pan)” 


- | (213) 


-B4 


ее ea 
а 3 By, | 14721. 
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APPENDIX Н PA 


REYNOLDS NUMBERS AND TRANSITION 


A Reynolds number characterizing the boundary-layer development in 
regions 2 апа 3 of the shock tube is now developed. 


u 
us ч-т 


Та the notation of reference 2, a Reynolds number characterizing 


Region 


т, 


х (ug -u,)^ 
the boundary layer behind а shock is defined therein as Ве = S е 
e 


W 
Using the transformations indicated by equations (21) and arbitrarily 
basing v on free-stream conditions, the Reynolds number for region 2 
may be written as 


R Ж uo(u,c М Е) up/ ug (m) 
7 v2 gu uj Ug) 


where Е,т are the coordinates of a point in region 2. Consider Е,т 
to be a point on the characteristic line through x,t (preceding sketch). 
The Reynolds number for points along this line, as a function of Е, 16 
(using eq. (B2)) 


ug(x - Е) Up/Ug l +M - u,/ag 
MGM ын) ш 


ІҒ the transition Reynolds number is known, equation (H2) can be used to 
determine the values of Е at which transition occurs. The maximum 


F2 9778s n 
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Reynolds number in region 2 occurs at the contact surface (Е = Ey) and 


equals 
a4 x 12 2 u | 
ЕТИ 


Equation (H5) can be used as a general index as to whether the boundary 
layer (along the characteristic line) is primarily laminar or primarily 


turbulent. А plot of Бер ъ against М, for air in region 1 (with 


T4 = 520? В) is given in figure 6. 


The transition Reynolds number for the boundary layer behind a 
shock wave has not yet been established. In the absence of more accur- 
&te information, the transition Reynolds number for incompressible flow 


over a semi-infinite flat plate [Re = о(0.5х106) might be used to ев- 


timate the transition point behind a shock wave. That is, take 
Res = 0(0.5x108) ав а rough estimate for the transition Reynolds number. 


This is probably & conservative estimate for the strong-shock cases, 
since the large amount of heat transfer to the shock-tube walls may have 
a very stabilizing effect on the boundary layer. 


u 
Region 3 (2 Š <s s) 


For & point in region 5, the Reynolds number as defined in refer- 


Region 
5 


/ 
ж Бъзтъ 


х,5 


EX 


enee Z becomes (using eqs. (D2) and v based on the free stream) 


uz(a,v + E) uz/a 


(24) 
Va 1+ъ ал. 


Rez = 


=- =w А е то Nn ou GR i EE урм 
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For a point on the characteristic line influencing point x,t (preceding 
sketch), the Reynolds number becomes 


uzX uz/ag 1 + Ма + а 1 +M - а/а a 
3х  uz/24 5 + 245 E+ ( Mp - Us да (н) 


Уз 1+1 а, 1 + Mz х 1 + Mz + 843 11 


Rez = 
5 


At the contact surface (Е = Е), equation (Н5) has the value 


2 
Ты. за (=) с ЖаН Жа зе 1 | (H6) 


Ihe criterion Rez = 0(0.5x10°) might be taken as a rough estimate for 
the transition Reynolds number behind the expansion wave. 


The boundary-layer characteristics presented in appendixes D and Е 
indicate a discontinuity in the boundary-layer profile across the contact 
surface. (Note that Res v/ Rez x = уз/У2-) This discontinuity does not 

2 2 


actually occur, and the theoretical discontinuity thus represents a de- 
ficiency of the present method. The effect on shock attenuation is prob- 


ably small. 


6407 
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APPENDIX I 


GENERATION OF PRESSURE WAVES BY WALL SHEAR AND HEAT ADDITION 


The generation of pressure waves by body forces and by heat sources 
in one-dimensional flow is investigated, and the results are used to 
compare the method of the present report with that of reference 6. 


Wave Generation by Body Forces and Heat Sources 


Assume a uniform flow of pressure p, velocity u, and so forth, in 
a tube of constant cross section. This uniform flow is assumed to be 
slightly perturbed by weak body forces and heat sources. The equations 
of motion are (neglecting the possibility of mass sources, since this 
case is treated in eq. (1) 


^u Au) _ Ә/ур 
(За бад) = z а 
ZI CS 


ИИ 
sc PESE = 


where f(x,t) is the body force per unit volume acting in the 4x-direction 
and a(x,t) is the heat addition per unit volume per unit time. The зуш- 

bol As represents the entropy perturbation of a particle and is related 

to the pressure and density perturbations by 


As Ар „дю 
са о (I2) 


The perturbation at any point x,t can be shown to equal 


Ap _ Др V p^ 
p р р 
э L(a гв) 
а т\р P 
x (I3) 
XTE 
Cy pu Cot 


me my re ee ee gee eee Сч р -—-—-—-—————— -- - 
———————-—-—_- Horne cap тъ ep ne ——— —— — 6 
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where 
- Е x- Е 
Др" _ y а\&,® aku CL &--u at 
p 2ра(1 + M) Cpl 
x Е -X 
Ар” _ Y q 2,5 Е u f Е, % = a- u 
p 2ра(1 - M) ср а, ag 


The upper limit on the integral for Ap 15 +ә or -œ depending оп 
whether М <1 or М> 1. Since equations (1) and (IL) are linear, the 
solution for the case where mass sources are also present can be obtained 
by adding equation (2) to equation (15). Note that an elemental heat or 
mass source generates a symmetric pressure-wave pattern (positive pres- 
sure waves propagating in the downstream and upstream directions), while 
an elemental body force generates an antisymmetric pressure-wave pattern 
(positive pressure waves propagating in the downstream direction and neg- 
ative pressure waves propagating in the upstream direction). 


Comparison with Reference 6 


In effect, reference 6 uses equations (13) to find the shock atten- 
vation in a shock tube. The value of а is obtained by averaging, 
across the tube cross section, the heat transfer at the wall and the 
viscous dissipation in the boundary layer. The value of Е is found by 
averaging the wall shear across the tube cross section. Thus, if ау, 
is the heat transferred into the boundary layer per unit wall area, H, - 
is the net dissipation in the boundary layer per unit wall area, and т ; 
is the shear per unit агеа exerted by the wall on the fluid (taken to 
be positive in the 4x-direction), then 


а= р („+ н) -а (+) 
(т4) 


zd. = 
f= Ty dw 


>| ез 


Substituting equation (I4) into equation (15) and integrating along 
appropriate characteristic lines should yield the same results for shock 
attenuation as those which were obtained in reference 6 from a somewhat 
different viewpoint. However, these results are not in agreement with 
the results obtained from equation (4) (the latter being the basis of 
the present report). 


6<0? 
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Equation (4) is based on v while equation (I4) contains terms re- 
lating to the heat transfer at the wall, viscous dissipation, and wall. 
shear. But, from boundary-layer theory, it can be shown that v is 
dependent on terms related to the heat transfer at the wall, viscous 
dissipation, and the velocity profile (see, e.g., eqs. (рза) and (D3c)). 
The heat-transfer and dissipation terms in equations (14) and (4) play 
the same qualitative roles, and therefore, will not be discussed further. 
However, the wall-shear term in equation (Т4) in some cases has a sign 
opposite to that of the velocity-profile term in equation (4). Thus, 
these terms are not always in qualitative agreement, and this leads to 
discrepancies between the results obtained from equations (14) and (4). 


Consider, for example, the boundary layer in a shock tube for the 
case of negligible heat transfer and dissipation. For this case, у de- 
pends only on the velocity-profile term. Similarly, only the V term 
is retained in equation (4). The signs of v and T and of the re- 


sulting pressure perturbations in regions 2 and 5 of the shock tube are 
summarized in the following table: | 


Perturbations Perturbations 
based on v based on T 


(eq (4)) (eq (14)) " 


Thus, when heat transfer and dissipation аге neglected, а perturbation 
solution based on v differs from & perturbation solution based on Ты 


in regard to the sign of the Ар” waves in region 2 and the Apt waves 
in region 3. 
Tt may be concluded that attenuation solutions based on equation 

14) differ from those based on equation (4) in the following respects: 
1) The integration of equation (14) along characteristic line ab tends 
to overestimate its contribution to shock attenuation, and (2) the inte- 
gration of equations (14) along characteristic line be tends to under- 
estimate its contribution to shock attenuation. For the limiting case 
M, > 1, the integrations along ab and be are negligible, and equations 


(4) and (I4) give the same results for shock attenuation. 
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When the boundary layer is thin, it is obvious that the attenuation 
calculation should be based on equation (4). If the viscous shear af- 
fects the entire cross section (as for a long shock tube), there is no 
longer a core of potential flow and the solution should probably be based 
on equations (14) (or on the characteristics method of ref. 6 when non- 
linearities become important). The quantities т.» Чу; and H, should 


then be based on unsteady pipe flow rather than on thin-boundary-layer 
theory. Іп practice, the boundary layer is generally sufficiently thin 
to permit the use of equation (4). 
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(b) Wave diagram for perfect fluid. 


27277746 РАДИ ДРЕ VLL 
77 777? 


IA 
ШИ 


(с) Flow in shock tube with — fluid. 
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Figure 1. - Shock-tube pbenomena. 
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640% 


Figure 2. - Characteristic lines appropriate for study of attenuation in 
shock tube. 
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(а) Laminar case (eq. (15)). 


Figure 4. - Boundary-layer thickness at point b in region 2. (ав in region 1 1s віг at 520° В; wall ів assumed a per- 
fect conductor (Том = Ту). 
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(b) Turbulent case (eq. (#5)). 


Figure 4. - Concluded. Bourndary-layer thickness at point b in region 2. Gas in region 1 ta air at 520° В; wall is 
assumed a perfect conductor (т, чет). 
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Figure 5. - Attenuation in sir-air shock tube. T, 2 Те 520° В; Prandtl number а, 0.70; ratio of specific heats тү, 1..4. 
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Figure 6. - Concluded. 
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Theory, turbulent boundary layer 
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——— Theory, laminar boundary layer 
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—-—— Theory, turbulent boundary layer 
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(a) Pressure ratio 24/54; 4.061; shock Mach number relative to wall 
Mas 1.544. 


Figure 7. - Pressure behind shock weve as function of distance from dia- 
phregm. Air-air shock tube. Ту = T, = 520? В; hydraulic diameter 
d, 1/7. 
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tH — ——Theory, laminar boundary layer 
(eq. (6) with Тр = Ts y = T4) 
—-—- Theory, turbulent boundary layer 
(ref. 6, Cp = 0.0581 Ве-1/5) 
-------Тһеогу, laminar boundary layer 
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(c) Pressure ratio p "nmn ‚ 7.455; shock Mach number relative to wall 
М, 1.518. 


Figure 7. - Continued. Pressure behind shock wave as function of dis- 
tance from diaphragm. Air-air shock tube. Ту = Ty = 520° R; hydraulic 
diameter 4, 1/7. 
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НЕН Theory, turbulent boundary layer 
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(а) Pressure ratio P4/P1 › 17.915; shock Mach number relative to wall 

‚ - Mg, 1.792. | 

Figure 7. - Concluded. Pressure behind shock wave as function of dis- 
tance from diaphragm. Аіг-віг shock tube. Тү = Ty = 5200 В; hydraulic 
diameter а, 1/7. | a 
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tal conditions of reference 6 (air in region 1; turbulent boundary layer; 
TQ = To, y = 520° Ry Bh, 1/8 ft). 
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